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L'M 


Area  of  Rt^1  Fin 

Constant  used  in  formula  for  viscous  lift  on  hulls 

Projected  Area  of  Mth  Hull  Section 

Frequency  Dependent  Coefficient 

•  ■ 

Frequency  Dependent  Coefficient 

Frequency  Dependent  Coefficient 

Frequency  Dependent  Coefficient 

Frequency  Dependent  Coefficient 

Frequency  Dependent  Coefficient 

Frequency  Dependent  Coefficient 

Frequency  Dependent  Coefficient 

Cross  Flow  Drag  Coefficient  for  M^  Hull  Section 

Cross  Flow  Drag  Coefficient  for  Rth  Fin 

Lift  Coefficient  for  Rth  Fin 

Frequency  Independent  Coefficient 

Frequency  Independent  Coefficient 

Frequency  Independent  Coefficient 

Frequency  Independent  Coefficient 

Cross  Flow  Drag  of  Mth  Hull  Section 

Cross  Flow  Drag  of  R^  Fin 

Complex  Amplitude  of  Wave  Exciting  Force  in  Heave 
Total  Force  Vector  in  Frequency  Domain 
Force  Input  to  System  Transfer  Matrix 


*Note  -  Much  of  this  notation  is  taken  from  Reference  1  and  is  used 
here  without  change 
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System  Transfer  Matrix 
Frequency  Response  Function 

Frequency  Response  Function 

Frequency  Response  Function 
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5max 

6(KT) 


Constant  Mean  Forward  Velocity  of  Ship 

Local  Total  Velocity  at  the  Hull  Section 

Local  Total  Velocity  at  the  R^  Fin 

Wind  Speed  in  Knots 
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System 
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Input 
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ABSTRACT 


The  development  of  a  nonlinear  frequency  domain 
simulation  which  accomodates  frequency  independent 
nonlinearities  is  described.  Although  developed  around 
the  dynamics  of  a  particular  naval  craft,  the  simulation 
can  be  applied  to  a  general  situation  and  need  not  be 
restricted  solely  to  the  prediction  of  ship  motions. 

The  workability  and  usefulness  of  the  simulation  is 
demonstrated  by  simulating  three  scenarios  with  the 
SWATH  6A  craft.  Two  scenarios  are  open  loop;  the 
third  includes  automatic  control.  These  scenarios 
are  also  simulated  in  the  time  domain  using  a  general¬ 
ized  nonlinear  time  domain  simulation  which  is  based  on 
the  same  modeling  of  the  craft  dynamics  as  that  which 
is  used  for  the  nonlinear  frequency  domain  simulation. 

The  results  obtained  with  these  simulations  agree  for 
each  of  the  scenarios  simulated. 

When  the  results  obtained  with  the  nonlinear 
frequency  domain  simulation  are  contrasted  with  those 
obtained  with  a  linear  frequency  domain  simulation 
which  employs  equivalent  linearization  for  the  handling 
of  nonlinearities,  a  lack  of  agreement  becomes  noticeable. 

This  highlights  a  limited  usefulness  of  equivalent 
linearization. 

ADMINISTRATIVE  INFORMATION 

This  work  has  been  jointly  funded  by  the  General  Hydromechanics 
Research  Program  Task  Area  SR0230101,  element  number  61153N,  work  unit 
1572-145,  and  by  the  Ships,  Subs,  and  Boats  Program,  Task  Area 
ZF43-421.  The  latter  funding  was  authorized  by  the  Naval  Material  Command 
(08T)  and  administered  by  the  Ship  Performance  Department  High  Performance 
Vehicles  Program  (1507). 


INTRODUCTION 

Frequency  domain  analysis  is  extremely  useful  in  the  study  of 
dynamic  systems.  When  it  is  known  that  a  system  is  dominated  by  linear 
dynamics,  frequency  domain  techniques  can  be  applied  in  a  rather  casual 
and  routine  manner.  This  situation  changes  when  the  system  exhibits  non- 
linearities  which  cannot  be  ignored.  Very  often,  when  nonlinearities  must 


be  recognized,  some  form  of  linearization  is  attempted  or  the  frequency 
domain  approach  is  abandoned  altogether.  This  can  be  unfortunate  since 
linearization  is  not  always  carried  out  appropriately  (in  which  case  mis¬ 
leading  results  are  obtained),  and  abandoning  the  frequency  domain  approach 
means  relinquishing  a  very  powerful  tool. 

This  report  describes  a  technique  which  does  not  force  a  choice  of 
either  alternative.  The  technique  is  an  iterative  process  in  the  frequency 
domain  which  will  accommodate  single  valued  nonlinearities.  It  was 
developed  around  the  dynamics  of  the  SWATH  6A,  a  Smal 1-Waterplane-Area- 
Twin-Hull  craft.  However,  it  can  be  applied  to  any  situation  in  which 
nonlinearities  are  significant  and  must  be  accounted  for  in  the  analysis. 

COORDINATE  SYSTEM 

Only  heave  and  pitch  are  treated  in  this  report.  Nothing  is  lost  by 
this  restriction  except  unnecessary  complication  since,  in  the  conventional 
modeling,  the  combination,  sway-yaw-roll,  is  assumed  to  be  uncoupled  from 
the  combination,  pitch-heave,  and  further,  surge  is  assumed  to  be  un¬ 
coupled  from  any  other  degree  of  freedom*. 

Motions  and  forces  are  referenced  to  a  frame  that  is  a  right-handed 
orthogonal  Cartesian  coordinate  system.  This  frame  does  not  rotate  but 
translates  with  constant  velocity,  Uo,  and,  thus,  is  an  inertial  frame. 

The  velocity,  U  ,  is  the  mean  forward  velocity  of  the  craft.  The 
o”igin  of  this  frame  lies  in  the  plane  of  the  mean  free  surface  and  the 
plane  of  the  mean  position  of  the  longitudinal  centerplarve  of  the  craft. 

*It  must  be  emphasized  that  the  technique  described  in  this  report  can  be 
applied  to  a  six-degree-of-freedom  coupled  system.  The  modelino  chosen  for 
demonstrating  the  technique  just  happens,  by  convention,  to  be  a  rather 
highly  decoupled  system. 


The  X  0  Y  plane  coincides  with  the  plane  of  the  mean  free  surface;  the 
ooo  ~ 

XQ-axis  points  toward  the  bow  of  the  craft;  the  YQ-axis  points  toward  port; 
and  the  Z0-axis  points  up.  The  craft  makes  small  motions  (i.e.,  pertur¬ 
bations)  about  this  coordinate  system.  The  reference  point  on  the  craft 
is  that  point  which,  at  calm  water  equilibrium,  coincides  with  the  origin 
of  the  coordinate  system. 

SWATH  CRAFT  SYSTEM  DYNAMICS  MODELING 
A  system  modeling  that  is  employed  for  predicting  surface  ship  motions 
and  one  that  is  used  at  the  David  W.  Taylor  Naval  Sh'p  Research  and  De¬ 
velopment  Center  (DTNSRDC)  for  predicting  motions  of  SWATH  craft  consists 
of  a  set  of  equations  usually  referred  to  as  differential  equations  with 
frequency  dependent  coefficients.  Considering  the  vertical  plane  only, 
the  steady  state  equation  of  heave  is: 


<"+fl33>  VB33VC33VA35'B+B355+C358  ’  * 


The  equation  for  pitch  is: 


^I5+A55^9+B55®+C559+A53^o+B53Zo+C53Zo  =  Te  ^ 


(1) 


(2) 


In  the  use  of  equations  (1)  and  (2),  it  is  assumed  that  the  steady  state 
can  be  reached.  (M  +  A33).  B33,  A35>  i^,  (Ig  +  Agg),  B55,  A^  and  B^ 
are  functions  of  frequency.  C33,  C35>  C55  and  C53  are  frequency 


*  The  bar  over  a  symbol  indicates  a  complex  quantity.  Also  only  the  real 
part  of  the  forcing  function  is  assumed. 
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independent.  The  amplitudes ,F  and  T,  of  the  forcing  functions  are  complex 
to  allow  for  phase  relationships  with  respect  to  the  corresponding  wave. 
Since  the  equations  (1)  and  (2)  are  linear  steady  state  equations,  they 
generate  sinusoidal  responses  to  sinusoidal  inputs;  i.e.. 


z0<‘>  ■  zoae'jUt*  <3> 

e(t)  =  eae"ja)t  (4) 

Both  ZQa  and  ea  are  complex  to  allow  for  phase  relationships  with  respect 

to  the  forcing  functions.  Substitution  of  (3)  and  (4)  into  (1)  and  (2), 

gives  the  following: 

zoa  *  *  V“)T  <5> 

9a  *  H53(u)F  t  HjjUJT  (6) 

Here 


and 


33 

T35 


53 


55 


(<*>)  =  k4/ko 

(to)  =-K2/Kq 
(“) 

(«)  =  Kj/K0 


Kj(w)  =  -to2(M  +  +  ^33 

K2(w)  =  -to  A35  -  J<oB35  +  C35 

K3 (to)  =  "w  A53  '  J“B53  +  C53 
K4(w)  =  -w2(I5  +_A55)-jo)B55  +  C55 
_  £j(w)  K2(u) 

0  K3(tj)  K4(to) 


(7) 

(8) 

(9) 

(10) 

(ID 

(12) 

(13) 

(14) 

(15) 


*0nly  the  real  parts  of  the  responses  are  assumed. 
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For  a  particular  wave  frequency  both  F(w)  and  T(w),  the  amplitudes 
of  the  wave  induced  exciting  force  and  moment  associated  with  this  wave, 
are  known  and  the  responses  Z  (w)  and  e  (w)  are  easily  computed.  Of 

Oa  a 

course,  this  can  be  done  for  any  sum  of  a  finite  number  of  frequencies  by 
computing  the  response  to  each  frequency  and  summing  the  responses.  This 
is  what  is  meant  by  a  linear  system. 

FREQUENCY  INDEPENDENT  NONLINEARITIES 

A  next  level  of  complexity  when  considering  frequency  domain  analysis 
results  from  the  introduction  of  frequency  independent  single  valued  non- 
linearities.  If  a  linear  frequency  domain  analysis  is  to  be  used,  some 
type  of  linearization  must  be  carried  out.  Reference  1  includes  a  dis¬ 
cussion  of  the  application  of  equivalent  linearization  to  the  same  prob¬ 
lem  selected  for  demonstrating  the  nonlinear  frequency  domain  simulation 
described  in  this  report.  If  the  system  is  weakly  nonlinear  such  that  the 
output  associated  with  any  single  frequency  input  can  be  approximated 
by  a  sinusoid  of  the  same  frequency  and  if  its  amplitude  is  adjusted 
to  a  best  least  squares  fit  to  the  nonlinear  response,  the  application  of 
equivalent  linearization  is  valid.  However,  when  the  input,  as  is  often 
the  case  in  ship  design  problems,  is  a  seaway  represented  by  a  sum  of 
a  finite  number  of  sinusoids  of  differing  frequency,  and  the  nonlinearity 
involved  is  associated  with,  for  example,  the  square  of  the  sum  of  these 
sinusoids  the  use  of  equivalent  linearization  may  be  inappropriate. 
NONLINEARITIES  MODELED  IN  SWATH  6A  SIMULATIONS 

The  nonlinearities  that  are  incorporated  in  SWATH  ship  modeling. 
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as  given  by  Lee'1  ,  are  frequency  independent  and  single  valued.  They 
are  the  viscous  lift  and  drag  on  the  hull  sections  and  fins.  For  the 

J.  L. 

Min  hull  section  (see  Figure  1)  the  viscous  lift  and  drag  are  given  by 
LBM  =  2  pVBM  APBM  ao  slrtoBM^C0SaBM^  ^1€ 


1  2 

°BM  =  2  pVBM  APBM  CDBM  slnaBM I S1 naBM I 


Where 


iL.  2 

APBM  =  t^ie  Pr°Jec'tecl  area  ^  "  hull  section  (m  ) 

3 

p  =  the  mass  density  of  seawater  (1031.9  Kg/m  ) 

a„  =  0.07 
o 

X  L 

CnRM  =  the  crossflow  drag  coefficient  at  the  hull  section 
(=0.5  for  all  sections) 

VRM  =  the  total.relati ve  velocity  between  the  craft  and  the  water 
at  the  Min  hull  section  (m/sec) 

aBM  =  the  angle!  of  attack  of  the  flow  at  the  hull  section 

V  =  -U  cos  sine  -  zi7  (i)  sine+z  (i)  cose  (18) 

100  i  z0BM  i  *0BM 

V_  =  -U  s i  no-Z  cos  et-X.lDM6H-z<7  (i)  cose+z£v  (i)  sine  (19) 

2  0  0  °BM  i  Z0BM  i  X0BM 

VBM  ■  V1  +  U2  <20> 

aBM  =  tan"1  V2/V2  (21) 

Here  (i)  and  £  (i)  are  the  .ertical  and  longitudinal  components, 

r0BM  *0BM  th 

respectively,  of  the  wave  orbital  velocity  at  the  Min  hull  section  and 

X  L 

for  the  iin  wave  component.  XgBM  is  the  x-coordinate  of  the  location  of 

the  M  body  section.  UQ  is  the  constant  forward  speed  of  the  craft.  The 

viscous  lift  and  drag  at  the  Rtfl  fin  (there  are  four  fins)  are  given  by 
1  2 

LFR  =  2.pVFRAFR  CLFRa  SlnaFR  (22) 

°FR  =  \  pAFR  CDFR  VFR  S1naFR  ' S1  n<xFR '  *23^ 


♦References  listed  on  Page  27. 


Where 


V0  =  -U  cose+Z  sine-Ec,  (i)  sine+?£u  (i)  cose 
3  0  0  i  Z0FR  !  X0FR 

V.  =  -U  sin 0  -2  cose  +XnFD 6  +is7  (i)  cose+zcy  (i)  sine 


OFR  ^7  v 
UK  i  Z0FR 


i  a0FR 


v2  =  V2  +  V2  (26 

VFR  v3  4 

ctpR  =  tan_1V3/V4  (27! 

here  l7  (i)  and  L  (i)  are  the  vertical  and  longitudinal  components 
•°FR  *0FR  th 

respectively,  of  the  wave  orbital  velocity  at  the  Rtn  fin  and  for  the  i 

wave  component.  XgRR  is  the  x-coordinate  of  the  Rth  fin.  UQ  is  the 

constant  forward  speed  of  the  craft.  For  the  SWATH  6A 

Ci  cd  =  4-38  for  the  two  forward  fins  and  3.4  for  the  two  aft  fins 

a 

CQpR  =  1.2  for  all  fins 

2 

ApR  =  8.05  m  (for  R=l,  the  area  of  the  forward  port  fin) 

? 

=  8.05  m  (for  R=2,  the  area  of  the  forward  starboard  fin) 

2 

ApR  =  24.03  m  (for  R-3,  the  area  of  the  aft  starboard  fin) 

=  24.03  m2  (for  R=4,  the  area  of  the  aft  port  fin) 

In  addition  to  the  nonlinearities  identified  above,  there  are  non- 
linearities  associated  with  the  modeling  of  an  active  controller.  Figure 
2  shows  a  flowchart  for  the  implementation  of  the  controller  algorithm. 


The  fin  order  is  proportional  to  heave  velocity,  Zg(t),  where  the  constant 

of  proportionality,  y,  has  the  value  3.048°/m  sec-1.  The  fin  angle, 6,  is 

positive  for  trailing  edge  up.  The  order  is  clipped  at  6Mftx  =  15°  and 

passed  through  a  first  order  lag.  This  equation  is  implemented  using 
(41 

Tustin's'  ’  method.  The  explicit  implementation  is  shown  in  BLOCK  4  of 


Figure  2.  The  rate  at  which  the  fin  moves  is  limited  if  necessary,  to 

6^ax  =  10°/sec.  In  the  implementing  of  the  controller,  clipping  and 
rate  limiting  can,  of  course,  be  strong  nonlinearities. 
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Flowchart  for  the  Fin  Angle  Control 


JH  'II'IIMIW 


NONLINEAR  FREQUENCY  DOMAIN  SIMULATION 

The  overview  of  the  nonlinear  frequency  domain  simulation  given  below 
is  followed  by  a  detailed  description. 

What  is  required  is  a  solution  to  the  vector  equation 


Here 


X(«)  -  [h(J)1  (£eXT(“)+£hL(“)) 

(28) 

x(»)  =  i zo^)  e>)}T 

(29) 

|jH(J)]=  H35U) 

H53U)  fy-gU) 

(30) 

£ext(o))  =  {FU)ext  t(g))ext}t 

(31) 

Fnl(w)  =  (FU)nl  T(oj)nl}T 

(32) 

(w)  is  the  complex  frequency  spectrum  of  wave 

induced  hydrodynamic 

r,c.  c  i_EXT 

forces  and  moments,  which  are  independent  of  ship  motions,  combined  with 

all  other  applied  forces  and  moments.  £^L(to)  is  the  complex  frequency 

spectrum  of  all  frequency  independent  nonlinear  forces  and  moments.  Among 

these  are  the  lift  and  drag  on  the  hulls  and  fins  as  well  as  nonlinear 

forces  and  moments  associated  with  active  control.  F...  (w)  cannot  be 

-NL 

obtained  by  a  direct  frequency  domain  operation.  It  can,  however,  be 
obtained  from  a  time  domain  computation  of  the  nonlinear  forces  and  moments 
and  a  transformation,  via  an  FFT,  to  the  complex  frequency  domain. 
Initially,  the  craft  is  assumed  to  have  no  motion  except  for  a  constant 
forward  speed.  Having  applied  F^U)  onlylF^U)  is  set  to  zero  at  this 
point)  one  has  a  motion  vector  XU).  This  motion  vector  is  transformed 
via  FFT  1  to  its  time  domain  equivalent  to  obtain  X(t).  One  can  now  carry 

*The  under  bar  indicates  a  vector. 
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out  all  motion  dependent  nonlinear  force  and  moment  computations  including 
those  contributed  by  active  control.  The  forces  and  moments  obtained  are 

_  L. 

passed  through  an  FFT  and  one  obtains  F^(w),  say,  at  the  N  iteration. 

One  now  computes  the  motion  vector  at  the  (N+l)*^  iteration  as  follows: 

XU)N+1  =  IFUatF^U)  +  Fnl(o.)N)*  (33) 

The  iterative  process  is  continued  until  convergence  is  achieved.  In 

choosing  a  discrete  frequency  domain  representation,  the  frequency  spacing, 

Aw,  is  selected  to  insure  a  fundamental  period  (T  „  =  2tt/Aw)  in  the  time 

domain  large  enough  to  be  a  representative  portion  of  the  response  to  a 

stationary  random  input.  Also,  the  highest  frequency  represented  (w  = 

ma  x 

(N-I)Aw)  should  be  at  least  twice  the  highest  frequency  at  which  signifi¬ 
cant  response  is  expected. 

Figure  3  shows  a  flowchart  for  the  nonlinear  frequency  domain  simu¬ 
lation.  The  simulation  starts  with  an  initialization  block  (not  shown) 
in  which  a  sine  wave  approximation  to  the  seaway  and  the  wave  induced 
forces  and  moment  are  computed.  These  are  motion  independent  forces  and 
moments  which  are  combined  with  all  other  applied  forces  and  moments.  The 


sum  is  labelled  F£XT 


(w) . 


the  state  vector,  Xjt)  = 


Wave  orbital  velocities  are  also  computed  here, 
{  ZQ(t)  0  (t)}T  is  set  to  zero  £  0LD(w) 


is  set  to  zero  as  is  the  state  vector,  X^w).  Nonlinear  forces  and  moments 


which  can,  of  course,  include  linear  terms,  are  computed  in  BLOCK  I  in  the 


time  domain  for  one  cycle  of  the  lowest  non-zero  frequency  used  in  the  fre¬ 
quency  domain  representation  of  the  system.  BLOCK  II  transforms  the  nonlinear 


forces  and  moments  to  their  frequency  domain  description  via  a  Fast  Fourier 


Transform  (FFT).  These  are  added,  in  BLOCK  III,  to  F^U).  This  sum  gives  a 


*In  the  actual  implementation  of  the  algorithm  (W)N  is  modified  to  increase 
the  region  of  convergence.  -NL 
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BLOCK  I 


Figure  3  -  Flowchart  for  the  Nonlineat  Frequency  Domain  Simulation 
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frequency  domain  description  of  the  total  forces  and  moments*  F(w)  acting 
on  the  craft.  For  the  first  iteration,  i=o,  BLOCK  V  comes  into  play,  the 
matrices,  tufl  >  (juU  ,  and  [g(w£]  are  initialized  and  F^w),  is  set 
equal  to  the  recently  computed  £(<d)  and  saved.  In  BLOCK  VI,  £'(w),  the 
input  to  the  system  transfer  matrix,  [>(«)]  ,  is  computed.  This  compu¬ 
tation  is  the  actual  force  acting  on  the  system  only  when  convergence  has 
been  achieved.  The  state  vector,  Y(“)>  in  the  frequency  domain  is  computed 
in  BLOCK  VII.  The  simulation  has  two  loops,  an  inner  loop  over  the  index 
i  and  an  outer  loop  over  the  index  k.  If  during  any  iteration  the  conver¬ 
gence  criteria  in  BLOCK  VIII  are  met,  the  simulation  stops.  If  not, 

FqLD(oj)  is  set  equal  to  the  last  value  of  £(w).  £(w)  and  juX(w)  are  trans¬ 
formed  to  their  time  domain  description  via  an  FFT-1  in  BLOCK  X.  When  i 
equals  i  MAX,  the  k-loop  is  activated,  i  is  reset  to  zero,  k  is  set  to 
k+1,  and  X^w)  is  replaced  with  the  latest  values  of  X^).  This  has  the 
effect  of  restarting  the  iteration  process  with  the  old  values  of  X^U) 
and  F^U)  updated.  This  very  important  updating  was  found  to  speed  up 
convergence  significantly.  The  simulation  continues  as  before.  If  i=0 
BLOCK  XII  is  activated.  Computations  carried  out  in  BLOCK  XII  and  BLOCK 
XIII  lead  to  the  computation  in  BLOCK  VI,  of  £'(u).  Following  BLOCK  XIII, 
the  simulation  continues  as  before. 

SCENARIOS  SIMULATED 

Three  scenarios  were  simulated  to  demonstrate  applications  of  the 
nonlinear  frequency  domain  simulation.  In  the  first  scenario,  the  SWATH 
6A  is  proceeding  at,  nominally,  equilibrium  depth  and  at  a  speed  of  10 
knots  into  a  head  sea  state  6.  The  seaway  is  modeled  by  a  Pierson- 
floskowitz  energy  spectrum  which  in  turn  is  approximated  by  15  sine  waves 
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in  random  phase  to  one  another.  The  modeling  of  the  seaway  is  discussed 
in  Appendix  A.  A  second  scenario  is  the  same  as  the  first  except  an  active 
controller  has  been  added.  A  third  scenario  is  a  repeat  of  the  first  for 
zero  speed.  For  the  third  scenario  the  seaway  is  approximated  by  12  sine 
waves  in  random  phase  to  one  another. 

The  output  from  the  simulations  is  obtained  using  three  different 
appoaches.  One  approach  uses  the  technique  described  in  this  report; 
the  second  approach  uses  equivalent  linearization  and  a  linear  frequency 
domain  simulation;  and  the  third  approach  uses  a  generalized  time  domain 
simulation  which  does  not  require  a  linearization  of  the  nonlinearities. 

The  second  scenario  was  not  simulated  with  the  linear  frequency  domain 
simulation  since  that  simulation  will  not  accommodate  active  controllers. 

PRESENTATION  OF  RESULTS 

The  results  presented  here  and  obtained  using  the  nonlinear  frequency 
domain  simulation  technique  are  discussed  in  relation  to,  and  are  compared 
with,  results  obtained  with  a  Generalized  Nonlinear  Time  Domain  Motion 
Predictor  for  SWATH  Craft  (Reference  2)  and  with  the  Mono-Hull  or  Twin- 
Hull  Ship  Motion  Predictions  Computer  Program  (Reference  3).  In  the  former, 
the  full  power  of  the  nonlinearities  is  modeled  as  it  is  in  the  Nonlinear 
Frequency  Domain  Simulation,  and  consequently,  the  steady  state 
solutions  obtained  should  agree  with  solutions  obtained  with  the  Nonlinear 
Frequency  Domain  Simulation*.  The  results  obtained  with  the  program  of 
Reference  3  do  not  agree  since  this  program  treats  the  nonlinearities  via 
equivalent  linearization  and  one  cannot  expect  agreement. 

*0nly  steady  state  solutions  can  be  obtained  with  the  Nonlinear  Frequency 
Domain  Simulation. 


..  -aL-aJLL,:-. 
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Figures  4,  5,  6,  and  7  show  results  obtained  when  the  three  simula¬ 
tions  carry  out  scenarios  1  and  2.  Figure  4  shows  the  heave  response. 

Figure  5  the  heave  phase,  Figure  6  the  pitch  response,  and  Figure  7  the 
pitch  phase.  In  all  these  figures,  the  responses  obtained  with  the  Non¬ 
linear  Frequency  Domain  Simulation  agree  with  those  obtained  with  the 
Generalized  Nonlinear  Time  Domain  Motion  Predictor  while  the  results 
obtained  with  the  linearized  frequency  domain  simulation  do  not  agree 
with  the  other  two  over  the  full  frequency  span.  It  is  important  to 
observe  that  in  Figures  4  and  6  both  the  Nonlinear  Frequency  Domain 
Simulation  and  the  Generalized  Nonlinear  Time  Domain  Motion  Predictor  show 
motion  amplitudes  at  frequencies  lower  than  0.45  rad/sec,  the  lowest  fre¬ 
quency  wave  input.  Of  course,  the  linearized  frequency  domain  cannot  and 
does  not  predict  such  motions.  Since  the  phase  plots.  Figures  5  and  7,  show 
phases  relative  to  the  input  wave  there  can  be  no  entires  on  these  figures 
for  frequencies  lower  than  0.45  rad,  the  frequency  of  the  longest  wave  input. 

Figures  8,  9,  10,  and  11  show  results  obtained  when  the  three  simu¬ 
lations  carry  out  scenario  3.  Since  the  forward  speed  for  this  scenario 
is  zero,  no  results  with  a  controller  are  included.  At  zero  speed,  the 
fins  are  not  force  effectors.  Again,  the  results  obtained  with  the  non¬ 
linear  frequency  domain  simulation  and  the  generalized  time  domain  simu¬ 
lation  are  in  agreement  and,  again,  motion  is  predicted  at  frequencies 
lower  than  the  lowest  frequency  wave  input.  The  results  obtained  with 
the  linear  frequency  domain  simulation  are  not  in  agreement.  In  fact, 
here  the  disagreement  is  more  significant  than  that  of  scenario  1.  In 
scenario  1,  the  disagreement  at  heave  resonance  is  approximately  11% 
whereas  in  scenario  3  the  disagreement  is  approximately  35%. 
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HEAVE  PHASE  (DEGREES) 


•  RESULTS  OBTAINED  WITH  GENERALIZED  NONLINEAR  TIME  DOMAIN  MOTION 
PREDICTOR  FOR  SWATH  CRAFT  (NO  CONTROLLER)  (REFERENCE  2) 

O  RESULTS  OBTAINED  WITH  NONLINEAR  FREQUENCY  DOMAIN  SIMULATION 
(NO  CONTROLLER) 

A  RESULTS  OBTAINED  WITH  CURRENT  DTNSRDC  LINEAR  FREQUENCY  DOMAIN 
SIMULATION  (REFERENCE  4) 

0  RESULTS  OBTAINED  WITH  THE  TIME  DOMAIN  SIMULATION  INCLUDING  A 
CONTROLLER 

□  RESULTS  OBTAINED  WITH  NONLINEAR  FREQUENCY  DOMAIN  SIMULATION  PL 
INCLUDING  A  CONTROLLER 


SCENARIO  1, 


'SCENARIO  2 


NOTE:  THE  SEAWAY  INPUT  IS 
A  HEAD  SEASTATE  6  AS  REPRESENTED 
BY  A  PIERSON  MOSKOWITZ  ENERGY 
SPECTRUM.  (SEE  FIGURE) 
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Figure  5  -  Heave  Phase  for  Scenarios  1  and  2 
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•  RESULTS  OBTAINED  WITH  GENERALIZED  NONLINEAR  TIME  DOMAIN  MOTION 
PREDICTOR  FOR  SWATH  CRAFT  (NO  CONTROLLER!  (REFERENCE  21 
O  RESULTS  OBTAINED  WITH  NONLINEAR  FREQUENCY  DOMAIN  SIMULATION 
(NO  CONTROLLER) 

A  RESULTS  OBTAINED  WITH  CURRENT  DTNSRDC  LINEAR  FREQUENCY  DOMAIN 
SIMULATION  (REFERENCE  4) 

0  RESULTS  OBTAINED  WITH  THE  TIME  DOMAIN  SIMULATION  INCLUDING  A 
CONTROLLER 

□  RESULTS  OBTAINED  WITH  NONLINEAR  FREQUENCY  DOMAIN  SIMULATION 
INCLUDING  A  CONTROLLER 


Figure  6 
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PITCH  PHASE  (DEGREES) 


0 


•  RESULTS  OBTAINED  WITH  GENERALIZED  NONLINEAR  TIME  DOMAIN  MOTION 
PREDICTOR  FOR  SWATH  CRAFT  (NO  CONTROLLER)  (REFERENCE  2) 

O  RESULTS  OBTAINED  WITH  NONLINEAR  FREQUENCY  DOMAIN  SIMULATION 
(NO  CONTROLLER) 

A  RESULTS  OBTAINED  WITH  CURRENT  DTNSRDC  LINEAR  FREQUENCY  DOMAIN 
SIMULATION  (REFERENCE  4) 

V  RESULTS  OBTAINED  WITH  THE  TIME  DOMAIN  SIMULATION  INCLUDING  A 
CONTROLLER 

□  RESULTS  OBTAINED  WITH  NONLINEAR  FREQUENCY  DOMAIN  SIMULATION 
INCLUDING  A  CONTROLLER 


•100 


NOTE:  THE  SEAWAY  INPUT  IS 
A  HEAD  SEASTATE  6  AS  REPRESENTED 
BY  A  PIERSON  MOSKOWITZ  ENERGY 
SPECTRUM.  (SEE  FIGURE) 
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Figure  7  -  Pitch  Angle  phase  for  Scenarios  1  and  2 
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Heave  (M) 


RESULTS  OBTAINED  WITH  GENERALIZED  NONLINEAR  TIME  DOMAIN  MOTION 
PREDICTOR  FOR  SWATH  CRAFT  (NO  CONTROLLER)  (REFERENCE  2) 

O  RESULTS  OBTAINED  WITH  NONLINEAR  FREQUENCY  DOMAIN  SIMULATION 
(NO  CONTROLLER) 


u?  (RADIANS/SEC) 

Figure  8  -  Heave  Amplitude  for  Scenario  3 
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HEAVE  PHASE  (DEGREES) 


Heave  Phase  for  Scenario  3 


•  RESULTS  OBTAINED  WITH  GENERALIZED  NONLINEAR  TIME  DOMAIN  MOTION 
PREDICTOR  FOR  SWATH  CRAFT  (NO  CONTROLLER)  (REFERENCE  2) 

O  RESULTS  OBTAINED  WITH  NONLINEAR  FREQUENCY  DOMAIN  SIMULATION 
(NO  CONTROLLER) 

A  RESULTS  OBTAINED  WITH  CURRENT  DTNSRDC  LINEAR  FREQUENCY  DOMAIN 
SIMULATION  (REFERENCE  4) 


NOTE:  THE  SEAWAY  INPUT  IS 
A  HEAD  SEASTATE  6  AS  REPRESENTED 
BY  A  PIERSON  MOSKOWITZ  ENERGY 
SPECTRUM.  (SEE  FIGURE) 


Figure  12  shows  a  time  history  of  one  cycle  of  fin  order  and  fin 
displacement  under  automatic  control.  Figure  13  shows  a  time  history  of 
one  cycle  of  heave  displacement  with  and  without  automatic  control. 

CONCLUSIONS 

The  technique  developed  in  this  report  for  handling  frequency  inde¬ 
pendent  nonlinearities  in  frequency  domain  analysis  is  shown,  by  example, 
to  be  a  viable  technique.  The  question  as  to  the  range  of  nonlinear 
mathematics  throughout  which  the  technique  can  be  successfully  applied 
has  not  been  answered.  This  is  an  important  question.  However,  the 
answer  to  this  question  requires  further  effort. 
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SEAWAY  DESCRIPTION 


The  time  history  of  the  surface  elevation  of  a  seaway  is  usually 
approximated,  in  digital  simulations,  by  a  finite  sum  of  sine  waves.  In 
this  report,  the  Pierson-Moskovitz  representation  of  the  seaway  is  used. 


P(u>)  is  the  energy  ordinate  at  frequency  w  in  meters 


u  is  the  frequency  in  rad/sec 

Vw  is  the  wind  speed  in  knots  =  18.82  x  'My^/Z 

H^3  is  the  significant  waveheight  in  meters 

A  plot  of  this  spectrum  for  a  Sea  State  6  (significant  waveheight 


of  4.57  m  is  shown  in  Figure  14  for  a  ship  moving  at  10  knots  into  a  head 


sea  and  0  knots  into  a  head  sea. 


Twelve  sine  waves  were  chosen  to  approximate  the  seaway  for  the  0 
knot  case  while  15  sine  waves  were  used  in  the  encounter  frequency  domain 
for  a  ship  moving  into  a  head  sea  at  10  knots.  For  the  former,  the  12 
sine  waves  are: 

.232sin  (0.3825t  +  Qj)  ,884sin  0.7255t  +  ng)  .418sin(1.0615t  +  ng) 

. 722s in  (0.4675t  +  n2)  .695sin (0.8075t  +  flg)  .351s1n(1.1476t  +  f?1Q) 

1.006s in  (0.5525t  +  flg)  .610sin (0.8925t  +  n7)  .296sin(1.2325t  +  fin) 

1.006sin  '0.6375t  +  «4)  .503s1n(0.9775t  +  fig)  . 253si n ( 1 . 3175t  +  0J2) 

and  for  the  10  knot  case,  the  15  sine  waves  are 

.198s in  (0.4675t  +  Yl)  .753sin  (0.8925t  +  y6)  . 433si n ( 1 . 3175t  +  Yll) 

.497sin  (0.5525t  +  y2)  .683sin  (0.9775t  +  y7)  .387sin ;i.3025t  +  y12) 

.722sin  (0.6375t  +  y3)  .613sin  (1.0625t  +  Yg)  .351sin(i.4875t  +  y13) 


.811sin  (o.7225t  +  y4)  .546sin ( l. I475t  +  Yg)  . 317si n ( i . 5725t 
,805sin  (0.8075t  +  y5)  .485sin(i.2325t  +  YlQ)  .287sin(i.6515t 


The  forces  and  moments  corresponding  to  each  of  the  15  sine  waves 
(per  meter  wave  amplitude)  are,  for  the  10  knots  case 

FORCES  (N) 

l287078sin  (0.4675t+0.0378+?^)  381785sin (1. 1476t+0.00Q4+Yg) 

1096017sin  !0.5525t+0.0727+Y2)  417333sin (1.2325t-0.0383+y 1Q) 

919726sin  (0.6375t+0. 1087+y3)  478880sin(1.3175t-0.0725+Yn) 

766055sin (0. 7225t+0. 1454+y4)  554157sin(1.4025t-0.0998+y12) 

651706s  i  n (0 . 8075t+0 . 1789+y5 )  630444si n( 1 . 4875t-0 . 1237+y 13 ) 

550677sin (0 .8925t+0 . 0504+y  6 )  698264si n( 1 . 5725t-0 . 1554+y 14 ) 

428417s i n (0 . 9775t+0 . 0436-ny  ? )  750755s i n( 1 . 6575t-0 . 2645+y 15 ) 

382384s i n ( 1 . 0625t+0 . 0328+Yg ) 

MOMENTS  (N-M) 

5132690sin  (0. 4675t- 1. 7093+y 2 )  12246073sin(1.1475t+1.4622+Y9) 

4604980s in  (0.5525t-l. 7184+y2)  13394644sin(1.2325t+1.4484+Y10) 

3277861sin(0.6375t-l.  7OO3+Y3)  13418266sin(1.3175t+1.4367+yn) 

1460695s i n ( 1 . 7225t- 1 . 6242+y4 )  12064673s i n ( 1 . 4025t+l . 4260+y  n ) 

885429s i n ( 0 . 8075t+l . 7 104+y  5 )  9508579s i n ( 1 . 4875 1+1 . 4161+y l3) 

4195827s i n (0. 8925t+l. 5856+Yg)  5830657sin(1.5725t+1.4152H'y ^4) 

7348714sin  (0.9775t+1.5098*y7)  1408949sin(1.675t+1.8505+Yl5) 

10133452sin  (1.0625t+1.4806+Yg) 

For  the  0  knot  case,  the  forces  and  moments  are 

FORCES  (N) 

1372533sin (0. 3825t+1.2741+n1)  540120sin (0.8925t+0.0113+fi7) 

1089057sin (0.4675t+0.0508+n2)  586240sin (0.9775t-0.0394+fig) 

84 34 39s in (0.5525t+0. 1015+a3)  601871sin ( 1.0625t-0.0698+fi9) 

64723 Is in (0.6375t+0. 1485+fi4)  497231sin (1. 1475t-0.0827+«1Q) 


FORCE  (N) (CONT) 


521593sin  (0.7225t+0.1665+Og) 
495974sin (0.8075t+0. 1553+Og) 

MOMENTS 

2128444sin  (0.38251-1.7810+Oj) 
1310672sin (0.4675t-1.6365+O2) 
662579sin  (0.5525t+0.9815+O3) 
3423026sin (0. 63751+1. 2786+04) 
60740 10s i n ( 1 . 7225t+l . 3157+Og ) 
7705720sin  (0.8075t+1.4132ufig) 


254421sin (1.2325t-0.0487+«11) 
62103sin (1. 3175t+2.5599+012) 
(N-M) 

7937930sin (0.8925t+1.4880+O7) 
4691892sin (0.9775t+1.4271+og) 
1976727sin (1. 06251- 1.5676+0^) 
9502340sin  (1.1475t-1.6881+O10) 
14473240sin  (1.23251-1.7145+O.j) 
13974946sin (1.3175t-1.7153+o12) 


DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  DTNSRDC  REPORTS,  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH¬ 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES.  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY,  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 
BASIS. 


